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ABSTRACT 


The  performance  of  two  seismic  event  detectors,  the  Fisher 
detector  and  the  conventional  power  detector,  was  compared.  Long-period 
records  with  low  signal-to-noise  ratios,  recorded  at  the  Alaskan  Long 
Period  Array,  were  used  to  establish  performance  characteristics,  sen¬ 
sitivities,  response  to  off-azimuth  events,  and  response  to  signal  distortion. 
Results  confirmed  the  theoretical  prediction  that  there  is  no  important  dis¬ 
tinction  between  the  detectors. 
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SECTION  I 
INTRODUCTION 

This  report  studies  the  performance  of  two  seismic  signal  de¬ 
fectors;  the  Fisher  detector,  (Edwards.  Benno,  and  Creasey.  1,67).  and  the 
beam  power  detector.  These  detectors  are  designed  to  automatically  process 
seismic  data  and  alert  the  analyst  whenever  some  criterion  is  met  which  in- 
dicates  the  presence  of  a  signal.  The  Fisher  detector  responds  to  similar- 
t.te,  among  the  outputs  of  the  array  elements,  while  the  power  detector  re¬ 
sponds  to  sudden  changes  in  the  average  power  level  over  the  array. 

This  report  investigates  the  response  of  the  detectors  to  sig¬ 
nals  of  varying  signal-, o-„oise  ratio;  their  performance  characteristics  in 
terms  of  false  alarm  rat.  versus  detection  probability;  ,heir  response  to  off- 
ztmuth  signals,  and  their  response  to  irregularities  in  the  data  such  as 

glitches,  spikes,  and  amplitude  and  phase  distortion.  Both  theoretical  and 
experimental  results  are  presented. 


I 
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I 
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SECTION  II 

DESCRIPTION  OF  THE  FISHER  DETECTOR 


The  Fisher  detector  is  an  algorithm  which  responds  to  similar¬ 
ities  in  the  time -delayed  traces  from  an  array  of  identically  oriented  seismo¬ 
meters.  T ime -de-laying  the  traces  has  the  effect  of  forming  a  beam  in  the 
direction  specified  by  the  delay  times.  A  time  averaging  feature  is  included 
so  that  signal  similarity  must  persist  for  some  time  before  the  detector  re¬ 
sponds.  Thus  brief  accidental  similarity  due  to  noise  will  be  suppressed.  In¬ 
coherent  noise  is  dissimilar  from  one  site  to  another,  whereas  propagating 
signals  are  similar.  Therefore,  the  detector  will  respond  to  propagating  sig¬ 
nals  of  sufficiently  great  duration. 

The  Fisher  detector  takes  advantage  of  the  fact  that  the  mean 
square  of  a  series  of  numbers  is  the  same  as  the  squared  mean  only  when  the 
numbers  are  identically  i.he  same.  The  numbers  in  this  case  are  the  time 
shifted  outputs  of  M  seismometers.  If  y.(t)  =  Xj(t -7^)  where  is  the  output 

of  the  i-th  seismometer  at  the  time  t,  and  r.  is  time  shift  required  to  form  a 
particular  beam,  the  detector  calculates  the  quantity: 
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-  V  y2 

M  Z_j  yi 


(t) 


M  2 

k  X  ^  I 


(i) 


i=l 


i=l 


The  Fisher  detec  or  is  the  inverse  of  ( 1 ),  normalized  by  multi¬ 


plying  by  (M 


[,  M  i 


A  time  average  of  the  sums  formed  is  taken, 


before  the  detector  output  is  formed.  Then  the  detector  output  is: 
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where  the  bar  denotes  averaging  over  some  time  before  the  output  time  t  . 

If  the  seismometer  outputs  are  ail  the  same  the  Fisher  output  increases  with¬ 
out  limit.  For  real  signals  with  noise  the  seismometer  outputs  are  not  exact¬ 
ly  the  same,  and  the  Fisher  output  reaches  some  finite  value  as  shown  below 
for  a  few  special  cases. 

Equation  2  has  the  form 


S(t)  = 


(M-l)  B2 


M 


i=l 


where  B  is  the  ordinary  beam  output  i.  e.  , 

M 

B(t)  *  T?  L  yi<t) 

i=  1 


The  conventional  beam  power  detector  is: 


(3) 


(4) 


P(t)  =  —V  B2  (5) 

<T  L 
O 


where  <t  is  the  mean  square  value  of  the  noise  preceding  the  signal  at  a 
single  site. 


II  -  2 


To  find  the  sensitivity  of  the  Fisher  and  conventional  power  de¬ 
tectors  we  assume  that  the  signals  of  interest  are  long  period  and  are  thus 
'veil  dispersed.  Then  the  peak  amplitude  of  such  a  signal  v/ill  consist  of  very 
nearly  sinusoidal  motion.  We  therefore  consider  the  response  of  the  detec¬ 
tors  to  one  frequency  component  of  the  signals,  in  the  presence  of  noise. 


The  amplitude  of  the  signal  with  wavenumber  k  will  be  A  . 

The  signal  amplitude  will  be  assumed  constant  over  the  array  for  the  moment. 

Site  positions  are  denoted  by  X.  ,  and  there  are  M  sites.  The  noise  varies 

J 

from  site  to  site,  with  amplitude  n.  and  phase  <b  at  the  i-th  site.  Then 

J  J 

the  signal  vector  is: 
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The  noise  has  been  assumed  uncor  "elated  with  the  signal. 

The  beam  pow;r  is  obtained  from  the  data  matrix  by 

M 


B  =  LHTTHL  =  A  A*  + 


1  V 

m2  J  J 
j=1 


*  =  AA*  +  JL  or  2 
M  ° 


(8) 


where  Li(p)  is  the  look  direction  vector,  and  p  =  k  . 

i  p  •  X 


**»  *  i 
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(9) 


n-3 


The  beam  power  detector  output  is  then: 


AA*  +  -L  a1 
M  o 


1  2 

-  O' 
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(10) 


Inserting  these  results  in  equation  3  for  the  Fisher  detector 


(M 


S  = 


'  »  (AA*  +  -5f  <rc2) 


M  , 

T3  E(“  -A) 

j  =  l 


-  AA: 


1  2 

-  O' 

M  o 


(11) 


Cross  terms  between  the  signal  and  noise  in  the  sum  of  squares  average  to 
zero,  and  11  reduces  to 


S  = 


(M 


I  1  2  \ 

-  1)  (  AAV  +  M  °-n  ) 


o-  2  (  i  .  —  ) 

O  M  ' 


(12) 


Assuming  M>>1,  this  is 


1  +  M 


just  as  for  the  beam  power  detector. 


(13) 


To  find  the  response  of  the  beam  power  detector  to  an  off- 
azimuth  signal  we  make  p  different  from  k  in  equations  6  and  9.  The  re 
suit  is 
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i(k  -  p)  •  Xj 


_L  (T  ^ 

M  O 


The  first  term  is  the  array  response  and  is  a  function  of  the  array  coordinates. 
Its  maximum  occurs  when  k  =  p  . 

If  we  denote  the  quantity  called  the  array  response  in  14  by 


_  .  f  ^  i(  k  -  p  )  •  X 

“W-T  L  e 

M  l  ~ 

J  =  1 


and  rewrite  the  expression  11  for  the  Fisher  detector  in  the  presence  of 
an  oil-azimuth  signal,  we  find 

(M-1)(AA*  a(p)  +  ±  <ToZj 

S  =  - 5 - l -  (1 

A  A*  (1  -  a(p))  +  a-/  (1  .  -L) 


Both  the  Fisher  detector  and  the  beam  power  detector  have  their  maximum 
when  the  look  direction  coincides  with  the  signal  wave  vector.  It  will  be  shown 
below  that  as  the  look  direction  moves  from  the  signal  azimuth,  the  first  mini¬ 
mum  in  their  outputs  occurs  at  the  same  change  in  azimuth.  Therefore  their 
resolutions  are  the  same. 


p  is 


The  derivative  of  16  with  respect  to  the  look  direction  vector 

(M-l)AA*  {  AA*  (1  -a)  +  ^2(1  -  ~)  J  +  AA*  j  AA*  <v  +  -1 

- 

{AA*  (1  -«)  +  <rj-  (I  -  | 
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and  this  is  zero  only  when  da/dp  ,  the  derivative  of  the  array  response  with 
respect  to  look  direction,  is  zero.  But  the  beam  power  detector's  derivative 
with  respect  to  angle  is 


_d  p  _  AA  d  a 

d  p  d  p 

and  this  is  zero  at  the  same  point. 


(18) 


Next  we  consider  the  effect  of  irregularities  in  the  data.  Such 
irregularities  may  be  random,  such  as  glitches  or  spikes,  or  systematic, 
such  as  variations  in  response  from  one  seismometer  to  another. 

When  a  glitch  or  spike  is  present  we  can  no  longer  use  signal 
vectors  to  represent  the  data.  These  irregularities  may  be  thought  of  as 
large,  roughly  constant  deviations  from  the  mean  at  one  seismometer.  They 
differ  from  one  another  in  their  time  duration  -  spikes  are  short  and  glitches 
are  long.  In  either  case  the  seismometers  outputs  will  no  longer  be  similar 
to  one  another,  and  the  Fisher  detector's  output  should  go  to  that  expected 
for  uncorrelated  noise. 

We  will  represent  a  glitch  as  a  site  whose  amplitude  is  con¬ 
stant  .it  G  over  the  time  averaging  interval.  Using  4  for  the  beam  output 
we  have  in  this  case. 


B2(t)= 


_1_ 
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M-  1 
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j=i 


(19) 


Carrying  out  the  squaring  gives 
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Alien  the  time  average  is  taken  the  cross  term  between  the  glitch  and  the  rest 
of  the  seismometers  goes  to  zero. 

The  first  term  may  be  represented  (almost)  as  a  beam  formed 
from  M-l  seismometers,  differing  from  that  quantity  only  by  a  factor  of 

(M-l)/M.  For  large  M  this  will  make  no  important  difference.  Thus  the 
beam  power  output  is 


b2|,)  =  b2m-i<'>  +  \ 

M 


(21) 


where  ^  is  the  "reduced"  beam.  Thus  we  can  conclude  that  the  be 

power  detector  will  always  be  larger  in  the  presence  of  a  glitch  than  in  it 
absence,  whether  a  signal  is  present  or  not. 


am 

s 


Next  we  insert  these  results  into  equation  3  for  the  Fisher 


de  .ector 


S(t)  = 


,m-»  K-, +  5 

M 


2  , 


-  T  y2  - 

M  Yj 

j  =  l 


(22) 


B2  J_  i 

M-l  M  ”  M 


It  can  be  seen  that  this  output  is  almost  the  same  as  that  without  the  glitch, 

except  that  (M  -  l)  G  /M2  has  been  added  to  the  numerator,  and  G2  (1  -  J-)  to 

M  M 

the  denominator.  A  little  algebraic  manipulation  shows  that  this  always  has 
the  effect  of  reducing  the  output.  The  limit  for  G  very  large  is  1,  as 
expected. 

Spikes  may  be  treated  in  very  much  the  way  glitches  are,  ex¬ 
cept  that  since  they  are  of  short  duration  they  have  much  less  effect  on  the 
output. 
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Next  suppose  that  the  response  to  a  plane  wave  varies  from 
station  to  station,  whether  due  to  variations  in  siesmometers  or  geology. 
The  signal  vector  is 


i(<f).  +  k 


v,  he  re  the  A's  and  the  s  differ  from  site  to  site.  Then  the  beam  power  i- 


M  M 

i‘E  L 

j=i  k-i  J 


Defining  A  ^  A^  and  6.  -  A.  -  A  we  have 

j  =  l  J 


M  M 


i?  E  El 

•  i  i  i  ' 


AA  -  A  8.  +  A  8 
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j  =  l  k=l 


' .  +  6.6  i  e 

J  J  kj 


-V  I  2 

"KT  %  (25) 


Tlie  second  sum  averages  to  zero,  as  do  those  terms  in  the  third  for  which 
j  ?  k.  Defining 


M  M 


EE 

j=l  k=l 


i  i<t>.  -  4>.) 

J  k 


2  1  1 


P  AA 


h  E 


\h 


th<  beam  power  detector  becomes 


1  7^2 


A_A  (cH  £-) 
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The  parameter  n2  is  a  measure  of  phase  distortion  and  ts 


near  one  for  a  practical  array.  Amplitude  distortion  is  measured  by  ,,2 

_ 11  .  .  -2  7  b 


R°aS °nable  vaklcs  P  He  between  .  01  and  .  02  at  ALPA 


The  expression  for  the  Fisher  detector  is 


S  = 
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JM-1)  j  AA  (a2  +  )  +  -£  <r2  j 


k  1 


j  AA  (a2  +  £.)  4  _!_  A 

(  M  M  o 


(29) 


which  becomes  after  time  averaging 
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2  ,  r.  12 
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(30) 


,xx  <"2  +  £>  +^-2 


(T 


Averaging  over  sites  as  before  leads  to 

RZ 

2  P 


S  = 


(M-n  <"  +  WD  ^ 


C7 
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AA  (\.  aZ  4  ft2(\  -  -Li 

M ; 


+  (l  -  rr) 


(31) 
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hrorn  these  expressions  it  can  be  seen  that  the  beam  powc 
output  increases  tndefinately  as  the  siBnal  amplitude  increases,  but  that  the 
Fisher  detector  saturates  at  a  level  dependent  on  the  kind  and  decree 


of 


distortion  present.  For  pure  phase  distortion,  the  Fisher  detector  output 


an  go  no  higher  than 


limited  to  the  value 


i-<t 


2  •  W"  only  amplitude  distortion  is  present,  it  is 


(M-l)  (£-4  1) 


P2  (1  -  — ) 
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II- 9 


.w  tu.i....,-  - 


- - - 


-- - * - — - 


.  0,1. 


When  both  are  present,  saturation  occurs  at 
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SECTION  III 

ANALYSIS  OF  THE  DATA 


The  purpose  of  this  section  is  to  test  the  predictions  of  Section 
II  with  real  data,  and  to  evaluate  the  performance  of  both  detectors  in  terms 
of  an  objective  criterion,  the  performance  characteristic.  Data  used  weie 
long-period  signals  from  events  in  the  3.  7  <  mb  <  4.  1  range  and  long-period 
noise  recorded  at  the  Alaskan  Long  Period  Array  (ALPA).  Processing  used 
about  1  3  sites. 

A.  FALSE  ALARM  RATE,  PROBABILITY  OF  DETECTION,  AND  T  IE 
PERFORMANCE  CHARACTERISTIC 

1.  Definitions 

The  program  output  takes  the  form  shown  in  Figure  HI-1.  ^ 

Time,  at  32  second  intervals,  is  printed  on  the  left.  Next  comes  a  row  of 
symbols,  one  for  each  azimuth.  These  symbols  represent  the  output  level  of 
the  Fisher  detector  for  the  first  time  gate,  at  the  time  in  question.  The  time 
gate  is  the  interval  over  which  data  are  averaged  in  the  detector  definition. 

The  range  of  outputs  which  gives  rise  to  a  particular  symbol  is  displayed  at 
the  end  of  the  plot,  as  part  of  a  histogram  which  gives  the  number  of  times 
each  symbol  occurred  for  each  time  gate  and  azimuth.  Next  comes  the  con¬ 
ventional  power  detector  output  levels  for  the  same  azimuths,  output  time, 
and  time  gate.  The  second  row  gives  the  Fisher  and  conventional  detector 
outputs  for  the  same  azimuths,  but  for  the  next  time  gate,  and  so  on. 

It  is  anticipated  that  the  detector  will  run  continuously,  form¬ 
ing  beams  which  cover  all  azimuths  with  some  rather  small  angular  spacing. 
Whenever  the  output,  on  a  particular  beam,  exceeds  some  predetermined  test 
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Azimuths  Spaced  at  15°  Intervals  Around  Primary  Beam  Direction 
Detector  Level  as  in  Table  III- 2 


FIGURE  III- 1 

PROGRAM  OUTPUT  FOR  15  BEAM  DIRECTIONS 
AND  TWO  TIME  GATES 
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level,  the  data  are  further  analyzed  to  determine  whether  the  output  was  due 
to  a  seismic  event,  or  to  some  form  of  noise.  The  detector  will  therefore 
function  as  a  flag,  which  alerts  the  analyst  to  the  need  for  further  investiga- 

v 

tion. 

The  detector  output  cm  reach  a  given  test  level  by  chan.e, 
even  if  the  input  is  noise,  and  this  is  called  a  false  alarm.  The  number  of 
times  this  happens,  per  unit  time,  is  the  false  alarm  rate.  On  the  other  hand, 
it  may  not  reach  the  level  required  for  a  detection  even  with  a  signal  present, 
if  the  signal  is  small  or  distorted.  The  probability  of  reaching  a  given  test 
level  when  a  signal  is  present  is  called  the  probability  of  detection  for  that 
level.  The  probability  of  reaching  a  given  test  level  is  of  course  either  zero 
or  one  for  a  particular  event,  depending  on  its  signal-to-noise  ratio.  Never¬ 
theless,  an  average  detection  probability,  which  will  tell  us  how  likely  it  is 
that  some  as  yet  uninvestigated  event  will  trigger  the  detector  can  be  esti¬ 
mated. 

* 

From  a  plot  such  as  Figure  III- 1  for  pure  noise  the  false  alarm 
rate  for  each  detector  can  be  calculated  as  a  function  of  time  gate  length.  The 
false  alarm  rate  for  a  given  power  level  is  defined  as  the  average  number  of 
times,  in  a  standard  time  interval,  that  the  detector  output  increases  past 
that  level.  Any  number  of  consecutive  output  values  which  exceed  a  given 
level  is  counted  as  only  one  false  alarm  at  that  level. 

Thus  if  the  detector  output  increased  past  the  power  level  re¬ 
presented  by  J  in  Figure  III  —  1 ,  four  times  in  a  8192  second  noise  sample, 
the  false  alarm  rate  for  the  J-th  level  would  be  4  per  8192  seconds.  False 
alarm  rates  are  normalized  to  256  seconds  in  this  report,  so  here  the  false 
alarm  rate  would  be  4/(8192/256)  =  .  125  per  unit  time  interval. 

The  performance  characteristic  is  found  by  eliminating  the 
threshold  level  between  the  detection  probability  and  the  false  alarm  rate. 
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TABLE  III-] 
EVENTS  STUDIED 


m  =  4.  1 
b 


KUR-171  -22  AL 
K  AM-  032  -  1  0  AL 
KAM-059-12AL 
KUR-235-03AL 
KAM-063  00AL 

KUR-228-22AL 
KAM-  1  92  -  12  AL 
TIB-037-04  AL 
IRA-041  -16QC 
TIB-073- 18AL 

IRA-154-00/  \ 
TUR  -  1 60-  12  AL 
IRA-187-21  AL 
RYU-  1  92-00AL 
PAK-214-09AL 


n^b  =  3.9 


IRA-  186-09AL 
TUR  - 196-04  AL 
BIB  -2  32  -  04  AL 
KOM-Ol  1  -  08  AL 
KAM- 01  6-1  1  AL 

KOM-042-13AL 
IRA-029-09AL 
CHI-030-03AL 
CRS-059-  i5QC 
KAM  -  180-14AL 

KUR-153-00AL 
SIN-  1  54  -  06  AL 
TIB  -  1  95  -  05  AL 
MED-  1  99-03  AL 
TIB-202-10AL 


mb  =  3’7 

KAM- 033 -04  A  L 
IT  Y-036-03AL 
KUR-054-05QC 
KUR-049-  14QC 
LOM-054-10  AL 

WRS-056-22QC 
KUR-066-09AL 
KUR-070-06AL 
KOR-074  -  1  5  AL 
OKH-078-19AL 

MON-  1  53  -  1 1  AL 
TSI-  154-1 6AL 
KUR-162-19AL 
IIQ-168-23AL 
ERS-172-09AL 

TUR-175-04  QC 
YUG-  1 77  -  04  A  L 
HIN-  178-20AL 
CKZ-181  -00QC 
KAM-  1  97  - 1  3  AL 

CKZ-200-06AL 
KAM-206- 1 OAL 
R  YU-2  14  -05  AL 
KUR-231-23AL 
KUR-236-10  AL 

KUR-241  -09AL 
TIB-242  -23  AL 
MON- 24  4  -  1 7  AL 
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FALSE  DE^C’^A^f  FISHER  AN°  CONVENTION 
DETEC  .ORS  AT  32  SECOND  TIME  GATE 
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3.  Results  and  Discussion 

Performance  characteristics  for  the  Fisher  detector  were 
calculated  separately  for  events  of  =  4.  1,  3.  9,  and  3.  7.  To  determine 
if  the  performance  characteristics  differed  from  magnitude  to  magnitude, 
tests  on  the  equality  of  the  probability  of  a  detection  at  fixed  false  alarm 
rates  but  different  magnitudes  were  performed.  It  was  found,  that  the  dif¬ 
ference  observed  was  not  significant  at  the  75%  confidence  level.  The  rea¬ 
son  for  this  effect  is  that  events  of  the  same  bodywave  magnitude  have  differ¬ 
ing  signal-to-noise  ratios  at  ALPA  due  to  distance  and  radiation  pattern  ef¬ 
fects  on  the  signals  and  temporal  variations  of  the  noise.  For  this  reason, 

all  the  events  studied  were  lumped  together,  and  the  results  of  this  report 

are  valid  for  m  near  4.  0. 

b 

Figure  III-4  shows  the  performance  characteristics  for  the 
Fisher  detector,  for  time  gates  of  length  32,  64,  128,  and  256  seconds.  A 
test  on  the  equality  of  the  probability  of  detection  at  a  false  alarm  rate  of 
.  2  per  256  seconds  shows  that  the  observed  deviations  between  the  32  second 
gate  and  the  256  second  gate  are  real  at  the  85%  confidence  level. 

This  confidence  level  is  not  as  high  as  one  might  hope  for,  and 
might  be  increased  by  taking  a  larger  sample.  However,  the  regular  increase 
in  detection  probability  with  gate  length  at  a  given  false  alarm  rate  is  addi¬ 
tional  evidence  that  the  difference  is  real.  Thus  it  is  concluded  that  for  low 
signal-to-noise  ratio  events  the  detector  with  a  256  second  signal  gate  is  sup¬ 
erior  to  that  with  a  32  second  gate,  at  false  alarm  rates  of  .  2  per  256  seconds 
and  lower. 

raise  alarms  are  due  to  fluctuations  in  the  mean  value  of  Ihc 
noise,  and  the  standard  deviation  of  the  mean  power  of  uncorrelated  noise  is 
proportional  to  2/N,  where  N  is  the  sample  length.  1  hus  inc  reas  ing  the 
gate  length  decreases  the  false  alarm  rate.  On  the  other  hand,  signal-to- 
noi se  ratio  drops  off  slowly  with  time  gate  length  so  long  as  the  gate  length  is 
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less  than  the  signal  length.  Since  the  detector  responds  to  signal-to-noise 
ratio,  this  means  that  detection  probabilities  drop  off  slowly  with  time  gate 
length.  This  explains  the  observed  result,  and  suggests  that  still  longer 
gates  would  giv2  more  attractive  performance  characteristics,  up  to  a  point 
which  depends  on  signal  length  and  thus  on  magnitude. 

At  higher  false  alarm  rates  this  superiority  of  the  larger  gate 
disappears.  However,  because  practical  limitations  would  preclude  operat¬ 
ing  at  high  false  alarm  rates,  this  effect  is  largely  academic.  For  example, 
a  false  alarm  rate  of  .2/256  seconds  corresponds  to  about  3  per  hour,  which 
already  might  be  close  to  the  upper  limit  that  the  overall  detection  system 
could  tolerate. 

The  confidence  level  for  the  hypothesis  that  the  other  pairs  of 
the  four  gates  examined  are  different  is  too  low  to  claim  that  they  are  in  fact 
different.  However,  it  is  reasonable  to  assume  that  the  256  second  gate  works 
best,  and  attention  will  be  restricted  to  it  for  the  remainder  of  this  report. 

Next,  in  Figure  III  -  5 ,  the  conventional  and  Fisher  detectors 
performance  characteristics  are  compated  at  a  time  gate  of  256  seconds.  The 
largest  deviation,  occurring  at  about  .  7  false  alarms  per  256  seconds,  is  not 
statistically  significant.  It  is  concluded  that  there  is  no  significant  difference 
in  the  performance  characteristics  of  the  Fisher  and  conventional  detectors 
for  256  second  time  gates. 

B.  RESPONSE  TO  OFF-AZIMUTH  SIGNALS 

To  find  the  response  of  off-azimuth  events  for  each  detector, 
the  detector  outputs  were  plotted  in  Figure  III-6  at  a  fixed  time  as  a  function 
of  azimuth  for  a  magnitude  4.  1  event.  The  side  lobes  of  the  response  are  ob¬ 
scured  by  the  noise  so  the  angle  at  which  the  detectors  reach  a  minimum  is 
not  visible.  Rejection  of  off-azimuth  signals  is  the  same  over  the  range 
studied,  as  predicted  in  Section  II. 
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Azimuth  in  Degrees  From  Primary  Beam  Direction 


FIGURE  III  -  6 

AZIMUTHAL  RESPONSE  OF  FISHER  AND  CONVENTIONAL 
DETECTORS  FOR  256  SECOND  TIME  GATE 


C.  RESPONSE  TO  DISTORTER  SIGNALS 

figure  III- 7  shows  the  amplitudes  of  1  isher  and  conventional 
detectors  as  a  function  of  time  for  RAM  -  168-09AL,  a  magnitude  4.  1  event. 

The  ou  put  of  the  Fisher  detector  falls  somewhat  below  that  of  the  conven¬ 
tional  r  owe r  detector. 


Amplitudes  for  this  event  were  so  low  that  amplitude  distortion 
could  not  be  measured  from  the  individual  site  time  traces.  Therefore  ^ 
was  assumed  to  be  .  015,  a  typical  value  at  ALPA.  The  results  of  Section  II 
were  used  to  calculate  a  value  of  a  of  .  73. 


The  parameter  a  corresponds  to  an  average  cosine  of  the  phase 
distortion  across  the  array,  and  is  thus  related  to  an  average  time  delay 
anomaly.  The  dominant  period  of  motion  for  this  event  was  20  seconds,  and 
this  leads  to  an  average  time  delay  anomaly  of  about  2  seconds.  We  might 
expect  delay  anomalies  of  half  this  amount  due  to  the  fact  that  individual  sites 
can  oe  time-aligned  only  within  a  sampling  period  of  2  seconds.  Variations 
of  the  signal  azimuth  with  time  and  genuine  departures  from  plane  wave  propa¬ 
gation  presumably  can  raise  the  time  delay  anomalies  to  2  seconds. 

Such  saturation  due  to  phase  distortion  is  no  detriment  to  the 
detector's  utility.  For  practical  arrays  where  o'  is  near  one,  saturation 
will  occur  at  a  level  well  above  that  at  which  signals  are  further  investigated, 
fhus  saturation  will  not  change  the  performance  characteristics  within  the 
range  of  reasonable  false  alarm  rates. 


D.  RESPONSE  TO  SIGNAL  IRREGULARITIES 

!•  igure  JII-8  shows  the  result  of  a  large  glitch  at  time  6  hours 
35  minutes  24  seconds.  The  conventional  detector  reaches  a  value  larger 
that  the  upper  bin  limit,  and  the  fisher  detector  goes  to  a  value  smaller  than 
the  smallest,  bin  limit.  This  is  in  agreement  with  tin c  results  of  Section  II, 
vdiere  it  was  predicted  that  a  glitch  would  drive  the  conventional  detector  to 
very  large  values  and  the  Fisher  detector  to  1. 
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Time  Fisher  Conventional 

06.30.04  -------- — 17 . 

06.30.36 - 1?  I 

06. 31  .OP - 77 - 1?  . 

06.31  .40  -7  7  717--7  7  7-77  7 - 7 - 1? 

06.3?. 12  -YYYYY-ZYYY-7YY7YY7YY- 7YYY77YY  1? 

06.32.44  - 7YY  7YY-7YYY77YY  12 

06.33.  16  - XVWXX  X7  YX  XX  7  ?XX  17 

06.33.48 - XWWYYY-7YYY — YY  1? 

06.34.20  - YXX  7 - l? 

06.34.62 - 000  JJJKPPFPPFPF  17 

06.36.24 - ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  12 

06.35.66  - ♦♦♦♦♦♦♦♦♦*•♦♦♦♦•♦  17 

06.36.28 - ♦  ♦♦♦♦♦♦♦•»■♦♦<•«*•♦  12 

06.37.00  - ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  1? 

06.37.32  - ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  17 

06. 3ft. 04 - ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  17 

06.38.36  - ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  12 

06. 39. Oft - ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  17 

06.39.40  - MNNRPRP  VVVDUVVV  1? 

06.40.12  - WVMVmiVTTTSRRRR  12 

06.40.44  - nr>DRr  rr  pr  rc  F  rcc  ^ 

06.41.16 - rrcono^Fnnr»nnno  1? 

06. 41.48  - acmnnnFnrmnnnn  12 

06.42.20  - ftrroODOFf'nnnnnn  17 

06.42.62  - rr.ronnnFnnnnnnn  12 

06.43.24  - rrnr>r>nnFnnnnrnn  1? 

06.43.66  - rrr  nonof  cponnon  12 

06.44.28  - r  nnFFF^OrOPRrcF  1? 

06.45.00  - Tim - 7-77  11 

06.45.37  -  12 

06.46.04  -  17 

06.46.36  -  17 

06.47.08  -  17 

06.47.40  -  1  7 

Azimuths  Spaced  at  15°  Intervals  About  Primary  Beam  Direction 
Detector  Level  as  in  Table  III  -  2 
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RESPONSE  OF  FISHER  AND  CONVENTIONAL 
DETECTORS  TO  GLITCH  IN  DATA 
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If  a  Slgnal  had  been  present  at  the  same  time  as  the  glitch, 
the  F isher  detector  would  have  missed  it.  The  beam  detector  will  give 

extra  false  alarms  due  to  glitches,  but  these  false  alarms  can  be  investigated 
for  the  presence  of  signals. 
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SECTION  IV 


CONC  LUSIONS 


A  theoretical  analysis  of  the  Fisher  detector  and  conventional 
detector  has  been  carried  out,  assuming  only  that  the  noise  is  not  correlated 
with  itself  or  the  signal.  This  analysis  shows  that  the  Fisher  and  conventional 
detectors  should  respond  in  the  same  way  to  signals  and  should  have  the  same 
response  to  off-azimuth  signals.  The  conventional  detector  is  not  sensitive 
to  phase  or  amplitude  distortion,  but  the  Fisher  detector  output  is  reduced 
by  phase  distortion,  although  not  by  amplitude  distortion. 

Experimental  results  confirm  these  predictions.  The  experi¬ 
mentally  determined  performance  characteristic,  an  objective  measure  of 
how  well  a  detector  performs,  is  the  same  for  the  Fisher  detector  and  con¬ 
ventional  detector,  within  the  limits  imposed  by  the  data  sample.  The  length 
of  time  over  which  the  data  are  averaged  has  a  significant  effect.  Longer 
times,  up  to  at  least  256  seconds,  give  better  results. 

The  rejection  of  off-azimuth  signals  is  the  same  for  the  det¬ 
ectors,  as  predicted.  Evidence  of  saturation  due  to  phase  distortion  in  the 
Fisher  detector  was  found,  and  was  consistent  with  the  predictions  of  Section 
II.  The  behavior  of  the  detectors  in  response  to  glitches  in  the  data  was  in¬ 
vestigated,  and  the  prediction  that  the  Fisher  detector  does  not  respond  to 
glitches  was  confirmed. 

One  of  the  difficulties  with  the  conventional  beam  power  detec- 
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tor  in  the  past  was  that  an  accurate  estimate  of  the  average  noise  power,  cr^  , 
was  hard  to  find.  It  was  anticipated  that  the  Fisher  detector  would  prove 
superior  to  the  conventional  detector  because  it  does  not  require  such  an 


estimate.  The  lack  of  distinction  between  the  detectors  found  in  this  study 
implicitly  suggests  that  the  method  used  to  estimate  a  is  a  good  one. 

A  number  of  suggestions  for  further  study  can  be  made.  The 
data  base  should  be  expanded,  both  by  investigating  more  events  in  the  mag¬ 
nitude  range  of  this  report,  and  by  taking  a  wider  range  of  magnitudes.  The 
performance  characteristics  may  be  distinguishable  if  this  is  done.  Longer 
time  gates  should  be  used,  as  suggested  by  the  results  of  Section  III. 


IV -2 


SECTION  V 
REFERENCES 

Alsup,  S.  A.,  and  E.  S.  Becker,  1973,  Simultaneous  Three-Component 

Broad  Band  Earth  Noise  Structure  at  Very  Long  Period  Experiment 
Stations,  Special  Report  No.  10,  Texas  Instruments  Incorporated. 

Barnard,  T.  E.  ,  1973,  Simulated  On-Line  Adaptive  Processing  Results  Us¬ 
ing  Alaska  Long  Period  Array  Data,  Special  Report  No.  2,  Texas 
Instruments  Incorporated. 

Edwards,  J.  P.  ,  III,  S.  A.  Benno,  and  G.  Creasey,  1967,  Evaluation  of  the 
CPO  Auxiliary  Processor,  CPO  Special  Report  No.  5,  Texas  Instru¬ 
ments  Incorporated. 

Strauss,  A.  C.  ,  1973,  Final  Evaluation  of  the  Detection  and  Discrimination 
Capability  of  the  Alaskan  Long  Period  Array,  Special  Report  No.  8, 
Texas  Instruments  Incorpoi  ated. 


V-l 


APPENDIX  A 

CALCULATION  OF  THE  NOISE  POWER 

A  question  not  addressed  in  Section  II  is  the  calculation  of 

2 

°i>  *  n°i®e  power.  Ideally  this  power  should  be  calculated  over  an  inter¬ 
val  which  is  known  to  be  free  of  signals.  On  the  other  hand,  the  noise  field 
changes  slowly  with  time,  so  a  single  value  for  the  noise  power  cannot  be 
used  for  all  time.  An  attempt  is  made  to  satisfy  the  requirements  of  con¬ 
tinuing  adaptation  and  exclusion  of  signals  in  the  present  detector.  Previous 
work  (Barnard,  1973),  did  not  meet  these  requirements  and  give  rise  to  un¬ 
acceptable  results. 

The  data,  before  being  input  to  either  detector,  is  quality 
checked  for  dead  or  unreasonably  large  channels.  After  the  beam  is  formed, 
in  a  particular  direction,  the  average  noise  power  for  that  beam  is  calculated, 
over  basic  time  gates  of  32  seconds.  The  logarithms  of  these  powers  are 
taken,  and  exponentially  smoothed  into  past  time.  That  is,  powers  at  pre¬ 
vious  times  get  exponentially  less  weighting  than  do  powers  at  recent  times. 

The  logarithm  is  taken  before  smoothing  because  it  has  been 
found  ( Alsup,  1973)  that  noise  power  distributions  are  skewed,  in  that  powers 
larger  than  the  mean  occur  more  often  than  half  the  time.  Taking  logarithms 
before  smoothing  reduces  the  influence  of  these  powers. 

The  smoothing  process  described  above  is  repeated  every  32 
seconds,  so  that  the  average  power  is  updated.  Wherever  a  noise  power  ex¬ 
ceeds  the  current  smoothed  average  by  6  dB  or  more,  it  is  assumed  that  a 
signal  is  present,  and  the  updating  process  is  stopped,  to  be  continued  after 
the  signal  has  passed.  In  this  way,  long  term  drifts  in  the  noise  power  are 
allowed,  but  signals  are  not  included  in  the  noise  power  estimate. 
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